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Computersystem und Verfahren zur Berechnung cines pharmakokinetischen 
Verhaltens einer chemischen Substanz in Insekten 

Die Erfindung betrifft ein Gdmputersystem zur Berechnung des pharmakokinetischen 
Verhaltens einer chemischen Substanz in Insekten auf der Basis eines physiolo- 
gischen Modells sowie ein entsprechendes Verfahren und Computerprogrammpro- 
dukt. 

Aus dem Stand der Technik sind zahlreiche physiologisch-basierte Modelle fur 
Saugetiere bekannt (z. B, Chamick et al., J. Pharmacokin. Biopharm. 23, 217 
(1995)). Femer ist aus dem Stand der Technik ein physiologisch-basiertes Modell ftir 
Rauperi bekannt (Greenwood et al., Pestic. Sci. 30, 97 (1990)). Hierbei wird das In- 
sekt durch Konipartimente beschrieben, die jeweils ein einzelnes Organ des Insekts 
reprasentieren. Die "Verschaltung" der einzelnen Kompartimente ergibt sich aus der 
bekannten Physiologic von Raupen. Wesentliche Parameter dieses Modells sind die 
Ratenkoeffizienten fur inter-kompartimentellen Massentransport, die die Geschwin- 
digkeit der Verteilung bestimmen, und die Organ-VerteilungskoefFizienten, die das 
Konzentrationsverhaltnis zwischen dem jeweiligen Organ und der Haemolymphe, die 
der Blutflussigkeit im SSugetier entspricht, im thermodynamischen Gleichgewicht 
angeben. 

Aus dem Stand der Technik ist es bekannt, solche kompartimentellen Modelle einzu- 
setzen, um ein experimentell ermitteltes pharmakokinetisches Profil einer Substanz 
retrospektiv durch Anpassung der Raten- und Verteilungskoeffizienten zu beschrei- 
ben (z. B. Lagadic et al., Pestic. Biochem. Physiol. 45, 105 (1993) & Pestic. Bio- 
chem. Physiol. 48, 173 (1994)). 

Der Erfindung Hegt demgegeniiber die Aufgabe zugrunde ein Computersystem zur 
Vorhersage des pharmakokinetischen Verhaltens einer chemischen Substanz in 



Le A 36 125 



-2- 

einem Insekt zu schaffen, sowie ein entsprechendes Verfahren und Computerpro- 
grammprodukt. 

Die der Erfindung zugrunde liegenden Aufgaben werden mit den Merkmalen der 
unabhangigen Patentanspriiche jeweils gelost. Bevorzugte-Ausfiihn^ der 
Erfindung sind in den abhangigen Patentanspriichen angegeben. 

Die Erfindung ermoglicht es, auf besonders effiziente Art und Weise eine Vorhersage 
des pharmakokinetischen Verhaltens einer chemischen Substanz in einem Insekt zu 
berechnen. Insbesondere ermoglicht es die Erfindung, die Aufiiahme, Verteilung und 
Ausscheidung von chemischen Substanzen in Insekten auf der Basis physikochemi- 
scher Parameter abzuschatzen. 

Hierzu wird ein physiologisch basiertes pharmakokinetisches Simulationsmodell 
eines Insekts zur Vorhersage von Konzentrations-Zeit-Verlaufen einer cheniischen 
Substanz in den Kompartimenten des Insekts verwendet. Das Simulationsmodell 
beinhaltet zumindest einen von der zu untersuchenden Substanz abhangigen Para- 
meter. Der Oder die Parameter des Simulationsmodells werden fur eine bestimmte 
Substanz basierend auf einer oder mehrerer der physikochemischen Eigenschaften 
der Substanz vorhergesagt. - 

Bei diesen physikochemischen Parametem handelt es sich beispielsweise um die 
Lipophilic der Substanz beschrieben durch den Verteilungskoeffizienten zwischen 
Wasser und phospholipidmembranen oder den Oktanol-Wasser- Verteilungskoeffi- 
zienten, das Molekulgewicht oder die LiosHchkeit. Die relevanten physikochemischen 
Parameter der Substanz konnen entweder aus einfachen Experimenten bestimmt oder 
mittels an sich bekannter Verfahren wie QSAR oder neuronalen Netzen direkt aus 
dem Deskriptor der chemischen Struktur der Substanz ermittelt werden. 

Im letzteren Fall ist es mit dem erfindungsgemaBen Verfahren sogar mogUch, virtu- 
elle, d. h. noch nicht synthetisierte Substanzen im Hinblick auf ihre Aufiiahme- und 
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Verteilungseigenschaften in Insekten zu bewerten. Weiterhin konnen aufgrund des 
etablierten Zusammenhangs zwischen physikochemischen und pharmakokinetischen 
Eigenschaften allgemeine Kriterien zur Optimierung von Insektizidenwirkstoffen 
abgeleitet werden. 

Nach einer bevorzugten Ausfuhnmgsform der Erfindiing wird als substanzabhangiger 
Parameter des Simulationsmodells der Ratenkoeffizient des inter-kompartimentellen 
Massentransports verwendet, welcher proportional zum Produkt aus der Permeabili- 
tat fur die Substanz und der effektive Oberflache der Kompartimente ist. Vorzugs- 
weise gibt es also eiiien substanzabhangigen Parameter flir jedes Kompartiment des 
Simulationsmodells, der die Permeabilitat des betreffenden Kompartiments fur die 
Substanz imd die effektive OberflSche des betreffenden Kompartiments beinhaltet. 

Besonders vorteilhafl ist, dass zur Ermittlung der Parameter fur das Simulationsmo- 
dell basierend auf einer physikochemischen Eigenschafl der Substanz keihe weiteren 
experimentellen Untersuchungen notwendig sind, sondern das der oder die Parameter 
basierend auf der physikochemischen Eigenschafl der Substanz ermittelt werden 
kann. Diese Ermittlung erfolgt auf der Grundlage einer Datenbasis, die fiir verschie- 
dene Testsubstanzen zuvor experimentell ermittelt worden ist. Die Datenbasis bein- 
haltet die fiir die Testsubstanzen experimentell ermittelten substanzabhangigen Pa- 
rameter des Simulationsmodells sowie die physikochemischen Eigenschaften der 
Testsubstanzen. Diese Datenbasis wird zur PrSdiktion des oder der substanzabhangi- 
gen Parameter fiir eine zu untersuchende Substanz verwendet. 

Nach einer bevorzugten Ausfuhrungsform der Erfindung wird aus der Datenbasis 
durch eine lineare Regression eine Berechnungsfunktion gewonnen. Beispielsweise 
handelt es sich bei der Berechnungsfunktion imi eine Funktion der Lipophilic und 
des Molekulargewichts. Zur Vorhersage eines Parameters fiir das Simulationsmodell 
fur eine zu untersuchende Substanz mussen also lediglich die Lipophilic und das 
Molekulargewicht der zu untersuchenden Substanz mit der Berechnimgsfunktion 
ausgewertet werden, um den Parameter zu erhalten. Mit dem vorhergesagten Para- 
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meter kann dann eine konkrete Simulation von Konzentrations-Zeit-Verlaufen bei 
der Aufhahme und Ausscheidung der Substanz in dem Insekt durchgefuhrt werden. 

Altemativ konnen anstelle der durch lineare Regression gewonnenen Berechnungs- 
funktion auch andere an sich^bekannte^Pradiktionsverfahren verwendet werdenr* ^ 

Von besonderem Vorteil ist, dass nachdem einmal eine Datenbasis fiir eine be- 
stimmte Anzahl von Testsubstanzen ermittelt worden ist, ftir die Simulation des 
phamiakokinetischen Verhaltens weiterer Substanzen keine zusatzlichen Experi- 
mente mehr erforderlich sind. Dies ermpglicht es mit einem hohen Durchsatz "Kan- 
didaten" fur potenzielle tnsektizide hinsichtlich deren phamiakokinetischen Verhal- 
ten zu beurteilen. Dadurch ist eine erhebliche Beschleunigung bei der Auffindung, 
Entwicklung und Optimierung von neuen Insektiziden moglich. 

Im Weiteren werden bevorzugte Ausfuhrungsformeri der Erfindung mit Bezugnahme 
auf die Zeichnungen naher erlautert. Es zeigen: 

ein Blockdiagramm einer Ausfuhrungsform eines erfindimgsgemaBen 
Computersystems, 

ein Flussdiagramm einer Ausfuhrungsform des erfindungsgemaBen 
Verfahrens, 

ein physiologisch basiertes Model! einer Raupe als Grundlage fur ein 
erfmdungsgemaBes Simulationsmodell mit. einem substanzabhangigen 
Parameter, 

ein Diagramm zur Darstellung der fiir das Simulationsmodell der 
Figur 3 gewonnenen Datenbasis zur . Durchfuhrung einer linearen 
Regression, lun eine Berechnungsfimktion zu ermitteln. 



Figur 1 



Figur 2 



Figur 3 



Figur 4 
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Die Figur I zeigt ein Computersystem 100. Dabei kann es sich um einen ublichen 
Personalcomputer (PC), eine Workstation oder um ein Client-Serversystem handeln. 

Das Computersystem 100 hat em Simulationsmodell 102. Hierbei haiidelt es sich um 
ein physiologisch basiertes pharmakokinetisches Model! eines Insekts. Hierzu bildet 
das Simulationsmodell 102 die Kompartimente des Insekts ab, und ermoglicht so die 
Vorhersage von Konzentrations-Zeit-Verlaufen einer Substanz in den Koniparti- 
menten des Insekts. 

Zu dem Simulationsnaodell 102 gehoren physiologische Parameter, 104 die nur von 
der Art des zu beschreibenden Insekts abhangen, soy/ie ein oder mehrere substanzab- 
hangige Parameter 106. Zur Durchfiihrung einer Simulationsberechnimg von Kon- 
zentrations-Zeit-Verlaufen ist also die Eingabe des Parameterwerts des substanzab- 
hangigen Parameters 106 fiir die auszuwertende Substanz erforderhch. Ein konkretes 
Beispiel fiir eine Ausfuhrungsform des Simulationsmodells 102 wird weiter unten 
mit Bezugnahme auf die Figuren 3 und 4 naher erlautert. 

Das Computersystem 100 hat femer eine Datenbank 108. Die Datenbank 108 dient 
zur Speicherung einer Datenbasis, die auf der Grundlage von experimentellen Unter- 
suchungen des pharmakokinetischen Verhaltens von Testsubstanzen in dem Insekt 
gewonnen worden sind. Fiir jede der zuvor experimentell ujitersuchten Testsubstan- 
zen beinhaltet die Datenbank 108 den oder die Parameterwerte, die fiir die betref- 
fende Testsubstanz experimentell ermittelt worden sind sowie zumindest eine physi- 
kochemisGhe Eigenschaft der betrefifenden Testsubstanz. 

Diese in der Datenbank 108 gespeicherte Datenbasis bildet die Grundlage fur die 
Vorhersage eines Parameterwerts fur eine neue zu untersuchende Substanz in dem 
Pradiktionsmodul 110. Beispielsweise wird durch ein lineares Regressionsyerfahren 
aus der in der Datenbank 108 gespeicherten Datenbasis eine Berechnungsvorschrift 
gewonnen, 4ie es ermoglicht, aus einer physikochemischen Eigenschaft einer zu 
untersuchenden Substanz, wie beispielsweise deren LipophiUe oder deren Moleku- 
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largewicht, den fiir die Durchfiihrung der Simulation mit dem Simulationsriiodell 102 
erforderlichen substanzabhangigen Parameterwert zu gewinnen. 

Das Computersystem 100 hat femer ein Eingabe-/Ausgabe-Modul 112 zur Eingabe 
der physikocheniischen Eigenschaft einer-zu untereuchenden Substanzv-Femer^ 
uber das Eingabe-/Ausgabe-Modul 112 die Ausgabe der simulierten Konzentrations- 
Zeit-Verlaufe. 

Das Eingabe-/Ausgabe-Modul 112 kann mit einer Datenbank gekoppelt sein, in der 
sich Deskriptoren von real existenten oder virtuellen Substanzen, wie z.B. potenziel- 
len Insektiziden, und deren physikochemische Eigenschaften befinden. Bei noch 
nicht synthetisierten, virtuellen Substanzen kann die fur die Eingabe in das Com- 
putersystem 100 erforderliche physikochemische Eigenschaft direkt aus dem De- 
skriptor der chemischen Striiktur der Testsubstanz mittels an sich bekannter Verfah- 
ren wie QSAR oder neuronalen Netzen ermittelt werden. Im Fall yon virtuellen Sub- 
stanzen kann anstelle der physikochemischen Eigenschaft auch der Deskriptor ilber 
das Eingabe-/ Ausgabe-Modul 112 eingegeben werden, Auf der Grundlage des 
Deskriptors wird dann in dem Computersystem 100 selbst die physikochemische 
Eigenschaft als EingangsgroBe fiir das Pradiktions-Modul 110 ermittelt. 

Die Simulationsergebnisse werden tiber das Eingabe-/Ausgabe-Modul 112 z,B- in die 
Datenbank eingegeben, sodass zu einem spateren Zeitpunkt eine Bewertung dei- 
Simulationsergebnisse erfolgen kann. 

Die Figur 2 zeigt ein entsprechendes Flussdiagramm. In dem Schritt 200 wird eine 
physikochemische Eigenschaft einer zu untersuchenden Substanz eingegeben. Bei 
der physikochemischen Eigenschaft handelt es sich beispielsweise um die Lipophilic 
Oder das Molekulargewicht der Substanz. 

In dem Schritt 202 werden die Parameterwerte der substanzabhangigen Parameter 
des Simulationsmodells fiir die Substanz auf der Grundlage einer Datenbasis basic- 
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rend auf der physikochemische Eigenschaft der zu beurteilenden Substanz berechnet. 
Diese Datenbasis beinhaltet zuvor experimentell ermittelte Parameterwerte von ver- 
schiedehen Testsubstanzen. £)iese Berechnung erfolgt in dem Schritt 202. 

In dem Schritt-204 werden die in dem Sehritt 202- berechneten Parameter in das 
Simulationsmodell eingegeben. Dort erfolgt darm die Berechnung von Konzentrati- 
ons-Zeit-Verlaufen der zu untersuchenden Substanz in den Kompartimenten des In- 
sekts. In dem Schritt 206 werden diese Konzentrations-Zeit-Verlaufe ausgegeben und 
konnen dann bewertet werden. 

Zur Ermittlung einer mSglichst aussagekraftigen Datenbasis fur die Pradiktion der 
Parameterwerte ist es vorteilhaft, wenn die zur experimentellen Ermittlung der 
Datenbasis verwendeten Testsubstanzen chemisch moglichst divers sind. Vorzugs- 
weise soUte eine neue zu untersuchende Substanz, fur die die Parameterwerte vorher- 
gesagt werden soUen, innerhalb des von den Testsubstanzen reprasentierten Ver- 
suchsraums liegen. 

Die Figur 3 zeigt beispielhaft ein physiologisches Modell fur eine Raupe. Dieses 
physiologische Modell basiert auf dem von Greenwood et al. (siehe oben) angegebe- 
nen Modell. • 

Zu dem Modell gehoren die folgenden acht Kompartimente der Raupe: 
Haemplymphe (hi) 300, Fett-Korper (fb) 302, Muskel (mu) 304, Oberflache der 
Kutikula (cs) 306, Kutikula (c) 308, Darmwand (gw) 310 und Danninhalt (gc) 312. 
Eine zu untersuchende Substanz wird entweder topisch fiber die Kutikula 308, oral 
fiber den Darm oder durch direkte Injektion in die Haemolymphe 300 verabrpicht. 

Die Haemolymphe 300 dient, ahnlich wie das Blut bei Saugetieren, als Haupttrans- 
portphase zwischen den verschiedenen Organen. 
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Stoffwechsel findet mit der Rate kc in der Kutikula 308, mit der Rate khl in der 
Haemolymphe 300 und mit der Rate kgw in der Darmwand 310 statt. Die Ausschei- 
dung findet mit der Rate kgc statt. 

Die Haemolymplie 300 wird als frei zirkulierende Flussigkeit modelliert, welche in 
Kontakt mit den durch die Kompartimente beschriebenen Raupenorganen steht. 
Transportvorgange zwischen den Kompartimenten erfolgen durch passive Diffusion 
mit Permeation uber die Membranen als ratenlimitierenden Schritt. Die Ratenkoeffi- 
zienten Xj^ des Massentransports zwischen den Kompartimenten werden durch das 
Permeabilitats-Oberflachenprodukt Fx Ax und das Volumen Vx der Organe bestimmt: 

= ^^^^ (x G {c,mu,fb,nc,gw}) 
^x 

Der Gleichgewichtszustand wird dabei nach einer Zeit tx » Vx/ (Fx Ax), erreicht. 
Das Verhaltnis der Konzentrationen in den peripheren Kompartimenten im Gleichge- 
wichtszustand wird dabei durch die Verteilungskoeffizienten bezogen auf die 
Haemolymphe Kx := Kx^ii^ bestimmt: 

C 

Kx=— ^ (x€{c,mu,fb,nc,gw}) (2) 

Chi 

Bel Kx handelt es sich also um den Verteilungskoeflfizienten einer Substanz zwischen 
Haemolymphe und dem Organ x im Gleichgewichtszustand. 

Zusatzlich miissen die Verteilungskoeffizienten zwischen der Oberflache der 
Kutikula 306 und der Kutikula (Kc/cs) sowie auch der Verteilungskoeffizient 
zwischen der Darmwand 310 und dem Darminhalt 312 (Kgc/gw) bekannt sein. 

Das von Greenwood et al. bekannte Modell mit acht Kompartimenten ist Txa: Beruck- 
sichtigung von 
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verschiedeneh Formen der Verabreichung (oral, topisch oder durch Injektion 
in die Haemolymphe), 

^ einen " Open Loop"-Massentransport:von^Nahrung uber den Emahrungskanal- 
und 

die Zeitabhangigkeit der Organvolumina wegen des Larvenwachstiims modi- 
fiziert worden. 

Basierend auf dieseni biophysikalischen Modell der Figur 3 lasst sich fur jedes Organ 
X der Raupe eine Massengleichgewichtsbeziehung durch eine Differenzialgleichung 
beschreiben. Solche Massengleichgewichtsbeziehungen sind fur die Organe der 
Raupe im Anhang angegeben (Formeln Al bis A8). 

Zur Berechnung eines pharmakokinetischen Profils mittels der Massengleichge- 
wichtsbeziehungen der Organe der Raupe mussen insgesamt 23 Parameterwerte be- 
. kannt sein. Hierbei handelt es sich um 

die Organvolumina Yx (acht Parameterwerte sowie zusatzliche Parameter- 
werte zur Beschreibung deren zeitlicher Veranderungen aufgrund des Larven- 
wachstums), 

die Verteilungskoeffizienten Kx der Subtanz zwischen Haemolymphe 300 und 
peripheren Kompartimenten (c, mu, fb, nc und gw), 

die Verteilungskoeffizienten zwischen der Oberflache der Kutikula und der 
Kutikula sowie zwischen der Darmwand und dem Darminhalt, 



die Stofifwechsebratenkonstanten in der Kutikula, der Haemolymphe xmd der 
Darmwand, sofem die zu untersuchende Substanz metabolisiert wird, imd 
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Ratenkoeffizienten 9ix = [PA]xA^x fur den Massentransport zwischen den 
Kompartimenten (c, mu, fb, nc und gw). 

Die Verteilungskoeffizienten sowie die Permeabilitats-Oberflachenpfodukte urid da- 
mit die Ratenkoeffizienten hangen von physiologischen Parametem des Insekts ab 
sowie von den physikochemischen Eigenschaften der zu untersuchenden Substanz. 
Femer sind auch die Stofiwechselratenkonstanten substanzspezifisch. 

Die Organvolumina konnen durch an sich bekannte Methoden experimentell eraiittelt 
werden. Die Verteilungskoeffizienten Kx der Substanz zwischen Haemolymphe und 
Organ X konnen wir folgt berechnet werden: 



ICj^ /HO' 

Kx — - ^xin-P ~ fwater^x ^fat ffat^ ^ pmieih fpmtein^x (3) 

^hl/HjO 



wobei 

Kfat = Verteilungskoeffizient der Substanz im Gleichgewicht zwischen Wasser 
und Fett (LipophiUe) • 

Kprotein = Vertcilungskoeffizient der Substanz im ^Gleichgewicht zwischen Wasser 
und Protein 

fwater.x = Volumeuanteil vou Wasscr in dcm Kompartiment X, 
ffat,x = Volumenanteil von Fett in dem Kompartiment X, 
fprotcin, x = Volumenanteil yon Protein in dem Kompartiment X. 
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Die MembranafFinitat (MA) oder altemativ der Oktanol/Wasser-Verteilungskoeffi- 
zient (Po/w) der Substanz kann als Schatzung des Verteilungskoeffizienten Kfat ver- 
wendet werden. Der Verteilungskoeffizient Kproiein kann beispielsweise aus der Bin- 
dungskonstanten von Humanserumalbumin (Kd"^"^ in [mmol]) und dem Molekular- 
gewicht des Proteins (65 kDa) ermittelt werden: - 



HSA = -4sa/65 ^ (4) 



Die Formel 3 kaiin durch Kombination der Fett- und Proteinanteile in einem einzigen 
organischen Anteil vereinfacht werden: 



Kx = — mit Kx/HjO = fwater,x + 0""fwater,x) K^fat (5) 

^hl/HjO 



Die Ratenkonstanten des Massentransports zwischen Kompartimenten werden 
experimentell fur die Testsubstanzen ermittelt, indem sie mit experimentellen phar- 
makokinetischen Daten abgeglichen und mit den physikochemischen Eigenschafteri 
15 der Testsubstanzen korreliert werden. 

Approximativ kann angenommen werden, dass der Permeabilitatskoeffizient P flir 
die Substanz fur alle Organe x identisch ist, d. h. Px = P fur alle Organe x, mit Aus- 
nahme der Kutikula. Femer ist P proportional zur Lipophilie (Kfat = MA oder 
20 - Kfat = Ko/w, d. h. dem Verteilungskoeffizienten OktanolAVasser) und dem Membran- 
diffiisionskoeffizienten Dmem der Substanz: 

P^K^D™™ (6) 



25 



Mittels einer Exponentialbeziehung wird die Abhangigkeit des Membrandififusions- 
koefi&zienten von dem Molekulargewicht (MW) der Substanz beschrieben: 
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Dmem ^ MW""-- (7) 
Die Formeln 6 iind 7 konnen miteinander kombiniert werden, so dass 

resultiert, wobei a eine Konstante ist. Mittels Gleichung 1 konnen die Ratenkon- 
stanten des Massentransports zwischen den Kompartimenten wie folgt ausgedriickt 
werden: . . \ 

X,=a^K,,MW-«- (9) 

. X 

Durch Logarithmierung erhalt man daraus 

Log(:^) = Log(a^)-S™emLog(MW) . (lO) 

Die Figur 4 zeigt eine entsprechende doppelt logarithmische Darstellung von A.x/Kfat 
uber dem Molekulargewicht MW mehrerer Testsubstanzen fiir verschiedene Kom- 
partimente. Ebenso zeigt die Figur 4 durch lineare Regressionen ermittelte Aus- 
gleichsgeraden durch die Messpunkte. 

Mittels Unearer Regression lasst sich also aus der experimentell ermittelten Daten- 
basis die Steigung Smem bestimmen sowie auch der Achsenabschnitt olAx/Vx- Damit 
hat man eine Berechnungsvorschrift erhalten, die es erlaubt, allein aufgrund der 
Lipophilic (Kfat) und des Molekulargewichts MW einer neuen, auch virtuellen Sub- 
stanz die GroBe zu berechnen. Aufgrund der Formeln 8 imd 9 erhalt man so P • Ax. 
Auf der Basis dieses substanzabhangigen Parameters sowie der Gleichgewichts-Ver- 
teilungskoefifizienten Kx bzw. Kc/cs und Kgc/gw aus Gl. (3) bzw. (5) lasst sich so das 



(8) 
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Gleichungssystem im Anhang z. B. numerisch losen, und man erhalt die Zeit-Kon- 
zentrations-Verlaufe der zu untersuchenden Substanz in den einzelnen Kpmparti- 
menten. 



LeA36 125 



-14- 

Anhang 

Gleichungssystem fur das pharmakbkinetische Simulationsmodell der Raupe 



. , Ve3^ = -[PA]e(Ce3--^) (Al) 

, - Ve ^ = [PA]e (Ces -- ^) + [PA], (Cy - keCc (A2) 

dC C 

run = -Z[PA]x (Cm -f^) - kh,Ch, (x € {c,mu.fb,nc.gw}) (A3) 

■ V„,%^ = [PAW(Chi-^) (A4) 



(fb) V*b^ = [PAk(C^-^) (A5) 



, , V„e^ = [PA]„e(Ch,-f^) (A6) 
(nc) dt 

. ^ Vg^^ = [PA]g^(Ci,,-^)-[PA]gw(Cg^--^)-VCgw (A7) 

(gc) ^gc^ = [PA]gw(Cgw--^) (A8). 

. gc/gw 

Cx: Konzentration der Substanz in dem Kompartiment x, 

Vx: Volumen des Kompartiments x; das Volumen kann nach einer experi- 

mentell ennittelten Funktion zeitlich veranderlich sein (V x(t)), . 

[PA]x: Permeabilitat-Oberfiachenprodukt von Kompartiment x, 

Kx/yi Verteilimgskoeffizient zwischen den Kompartimenten x und y, 

Kx = Kx^j: Verteilungskoeffizient zwischen Organ x imd Haemolymphe, 

kxi Ratenkonstante fiir den Stoffwechsel (x = c, hi, gw) 
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Patentanspruche 

1. Computersystem zur Berechnung eines pharmakokinetischen Verhaltens einer 
chemischen Siibstanz in Insekten mit: 

einem physiologisch basierten pharmakokinetischen Simulationsmo- 
del (102) eines Insekts zur Vorhersage von Konzentrations-Zeit-Ver- . 
laufen der chemischen Substanz in Kompartimenten des Insekts, wo- 
bei das Simulationsmodell zimiindest einen von der Substanz abhan- 
gigen Parameter aufweist, . . 

einem Vorhersagemodul (110) zur Vorhersage des zumindest einen 
Parameters basierend auf einer physikochemischen Eigenschaft der 
Substanz. 

2. Computersystem nach Anspruch 1, wobei es sich bei dem Parameter um das. 
Produkt aus der Permeabilitat fur die Substanz und der effektiven Oberflache 
der Kompartimente handelt. 

3. Computersystem nach Anspruch 1 oder 2, wobei es sjch bei dem Parameter 
imi den Ratenkoeffizienten : fur inter-kompartimentellen Massentransport 
(A,x = P AxA^x) handelt, wobei insbesondere das Volimien zumindest eines der 
Organe des Insekts eine Funktion der Zeit ist (Vx = Vx(t)). 

4. Computersystem nach Anspruch 1, 2 oder 3, wobei es sich bei dem Parameter 
um den Gleichgewichts-Verteilungskoeffizienten zwischen Organ und 
Haemolymphe (Kx, x € {c, mu, fb, nc imd gw}), zwischen Oberflache der 
Kutikula und Kutikula (Kc/cs), oder zwischen Darmwand und Darminhalt 
(Kgc/gw) handelt. 
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5. Computersystem nach einem der vorhergehenden Anspriiche 1 bis 4, wobei es 
sich bei der physikochemischen Eigenschaft nm den Verteilungskoeflizienten 
zwischen Wasser und Phospholipidmembranen, den Oktanol-Wasser-Vertei- 
lungskoeffizienten, das Molekulgewicht, die Loslichkeit, und/oder eine Kom- 
^bination dieser Parameter der Substanz handelt. — ^ ^ . ^ 

6. Computersystem nach einem der vorhergehenden Anspriiche 1 bis 5, mit 
einem QSAR-Modul oder einem neuronalen Netz zur Bestimmung der physi- 
kochemischen Eigenschaft aus einem Deskriptor der chemischen Struktur der 
Substanz. 

7. Computersystem nach einem der vorhergehenden Anspriiche 1 bis 6, wobei 
das Vorhersagemodul auf einer Datenbasis (108) basiert, die die physiko- 
chemischen Eigenschaften von Testsubstanzen und experimentell fur die 
Testsubstanzen ermittelten Parameter beinhaltet. 

8. Computersystem nach einem der vorhergehenden Anspriiche 1 bis 7, wobei 
das Vorhersagemodul eine Berechnungsfunktion zur Berechnung des Para- 
meteirs aus der Lipophilic und/oder dem Molekulargewicht der Substanz be- 
inhaltet. ^ 

9. Computersystem nach Anspruch 8, wobei die Berechnimgsfunktion auf einer 
linearen Regression experimentell ermittelter Parameterwerte basiert. 

10. Verfahren zur Berechnung eines pharmakokinetischen Verhaltens einer 
chemischen Substanz mit einem physiologisch basierten pharmakokinetischen 
Simulationsmodell eines Insekts zur Vorhersage von Konzentrations-Zeit- 
Veriaufen der chemischen Substanz in Kompartimenten des Insekts, wobei 
das Simulationsmodell zumindestens einen von der Substanz abhangigen 
Parameter aufweist, mit folgenden Schritten: 
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Eingabe einer physikochemischen Eigenschafl der Substanz in ein 
Vorhersagemodui zur Vorhersage des zumindest einen Parameters fiir 
die Substanz, 

5 - . - - Durchfuhning einer Simulation .mit Hilfe des Siinula^^^ 

Vorhersage von Konzentrations-Zeit-VerlSufen der chemischen Sub- 
stanz basierend aufdeni vorhergesagten zumindest einen Parameter. 

11. Verfahren nach Anspruch 10, wobei es sich bei dem Parameter urn das Pro- 
A 10 dukt aus der Permeabilitat fur die Substanz und der effektiven Oberflache der 

Kompartimente handelt. 

12. Verfahren nach Anspruch 10 oder 11, wobei es sich bei dem Parameter mn 
den Ratenkoeffizienten fur inter-kompartimentellen Massentransport 

15 ( x = PAxA^x) handelt, wobei insbesondere das Volumen zumindest eines 

der Qrgane des Insekts eine Funktion der Zeit ist (Vx = Vx(t)). 

13. Verfahren nach Anspruch 10, 11 oder 12, wobei es sich bei dem Parameter 
um den Gleichgewichts-Verteilvmgskoeffizienten zwischen Organ und 

20 Haemolymphe (Kx, x e {c, mu, fb, nc und gw}), zv^ischen Oberflache der 

a^JlK Kutikula und Kutikula (Kc/cs), oder zwischen Darmwaiid und Darminhalt 

(Kgc/gw) handelt. 

14. Verfahren nach einem der vorhergehenden Anspruche 10 bis 13, wobei es 
25 sich bei der physikochemischen Eigenschafl um den Verteilungskoeffizienten 

zwischen Wasser und Phospholipidmembranen, den Oktanol-Wasser- Vertei- 
lungskoeffizienten, das Molekiilgewicht, die Loslichkeit, und/oder eine 
Kombination dieser Parameter der Substanz handelt. 
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15. Verfahren nach einem der vorhergehenden Anspriiche 10 bis 14, wobei die 
physikochemische Eigenschaft der Substanz ifiittels eines QSAR oder mittels 
einem neuronalen Netz bestimmt wird. 



5_ 16. Verfahren nach einem der vorhergehenden Anspriiche 10 bis 15, wobei-die 
Vorhersage des zumindest einen Parameters basierend auf physikochemi^ 
schen Eigenschaften von Testsubstanzen xmd experimentell fur die Testsub- 
stanzen ermittelten Parameterwerten erfolgt. 



^ L \ 10 17. Verfahren nach einem der vorhergehenden Anspriiche 10 bis 16, wobei die 

Vorhersage des zumindest einen Parameters mit einer Berechnungsfunktion 
aus der Lipophilic und/oder dem Molekulargewicht der Substanz erfolgt. 

18. Verfahren nach einem der vorhergehenden Anspriiche 10 bis 17, wobei die 
15 Berechnungsfunktion auf einer linearen Regression experimentell ermittelter 

Parameterwerte basiert. 

19. Computerprogrammprodukt, insbesondere digitales Speichermedium, mit 
Programmmitteln zur Berechnung eines pharmakokinetischen Verhaltens ei- 

20 ner chemischen Substanz mit einem physiologisch basjerten pharmakokineti- 

schen Simulationsmodell eines Insekts zur Vorhersage von Konzentrations- 
Zeit-Verlaufen der chemischen Substanz in Kompartimenten des Insekts, wo- 
bei das Simulationsmodell zmnindestens einen von der Substanz abhSngigen 
Parameter aufweist, mit folgendeh Schritten: 



25 



Eingabe einer physikochemischen Eigenschaft der Substanz in ein 
Vorhersagemodul zur Vorhersage des zumindest einen Parameters fur 
die Substanz, 
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Durchfiihrung einer Simulation mit Hilfe des Simulationsmodells zur 
Vorhersage von Konzentrations-Zeit-Verlaufen der chemischen Sub- 
stanz basierend auf dem vorhergesagten zumindest einen Parameter, 

^20. - - Computerprogrammprodukt nach Anspruch45,-wpbei-^es sich^^ 

meter um das Produkt aus der Permeabilitat fiir die Substaxiz und der effekti- 
ven Oberflache der Kompartimente handelt. 

21. Computerprogrammprodukt nach Anspruch 19 oder 20, wobei es sich bei 
dem Parameter um den Ratenkoeffizienten fur inter-kompartimentellen 
Massentransport (A,x = PAxA^x) handelt, wobei insbesondere das Volumen 
zumindest eines der Organe des Insekts eine Funiktion der Zeit ist 
(Vi = Vx(t)). 

22. Computerprogrammprodukt nach Anspruch 19, 20 oder 21, wobei es sich. bei 
dem Parameter um den Gleichgewichts-Verteilungskoeffizienten zwischen 
Organ und Haemolymphe (Kx, x e {c, mu, fb, nc und gw}), zwischen Ober- 
flache der Kutikula xmd Kutikula (Kc/cs), oder zwischen Darmwand und 
Darminhalt (Kgc/gw) handelt 

23. Computerprogrammprodukt nach einem der vorhergehenden Anspriiche 19 
bis 22, wobei es sich bei der physikochemischen Eigenschafl um den Vertei- 
limgskoeflizienten zwischen Wasser und Phospholipidmembranen, den Okta- 
nol-Wasser-Verteilungskoeflfizienten, das Molekulgewicht, die Loslichkeit, 
und/oder eine Kombination dieser Parameter der Substanz handelt. 

24. Computerprogrammprodukt nach einem der vorhergehenden Anspriiche 19 
bis 23, mit einem QSAR-Modul, oder einem neuronalen Netz zur Bestim- 
mung der physikochemischen Eigenschaft aus einem Deskriptor der 
chemischen Struktur der Substanz. 
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25. Computerprogrammprodukt nach einem der yorhergehenden Anspriiche 19 
bis 24, wobei das Vorhersagemodul auf einer Datenbasis (108) basiert, die die 
physikochemischen Eigenschaften von Testsubstanzen und experimentell fiir 
die Testsubstanzen ermittelten Parameter beinhaltet. 

26. Computerprogrammprodukt nach einem der vorhergehenden Anspriiche 19 
bis 25, wobei das Vorhersagemodul eine Berechnungsfunktion zur Berech- 
nung des Parameters aus der Lipophilic und/oder dem Molekulargewicht der 
Substanz beinhaltet. 

27. Computerprogrammprodukt nach einem der vorhergehenden Anspriiche 19 
bis 26, wobei die Berechnungsfunktion auf einer linearen Regression experi- 
mentell ermittelter.Parameterwerte basiert. 
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Computersvstem und Verfahren zur Berechnung eines pharmakokinetischcn 
Verhaltens einer chemischen Substanz in Insekten 

- - Z u s a m m e n f a s s u n g . . ^ 

Die Erfindung betrifft ein Computersystem zur Berechnung des pharmakokinetischen 
Verhaltens einer chemischen Substanz in Insekten mit: 

einem physiologisch basierten pharmakokinetischen Simulationsmodel (102) 
eiries Insekts zur Vorhersage von . Konzentrations-Zeit-Verlaufen der 
chemischen Substanz in Kompartimenten des Insekts, wobei das Simula- 
tionsmodell zumindest einen von der Substanz abhangigen Parameter auf- 
weist, 

einem Vorhersagemodul (110) zur Vorhersage des zumindest einen Para- 
meters basierend auf einer physikochemischen Eigenschaft der Substanz. 



(Fig. 1) 
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Eingabe physiko-chemische Eigenschaft Test-^ubstanz 



^00 



Beiechming Parameterweite Test-Substanz aus Datenbasis experimentell zuvor 

untetsuchter Substanzen 


- - - 




Eiqgabe Paiaznet 


srwerte in Modell 



Aiisgabe ICoiizentrationsvetteihing Test<^ubstaiiz in Abhangigkeit von der Zeit 



Fig. 2 
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